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Abstract. In the present paper, the Tsallis statistics in the grand canonical ensemble was reconsidered in a 
general form. The thermodynamic properties of the nonrelativistic ideal gas of hadrons in the grand canon¬ 
ical ensemble was studied numerically and analytically in a finite volume and the thermodynamic limit. 

It was proved that the Tsallis statistics in the grand canonical ensemble satisfies the requirements of the 
equilibrium thermodynamics in the thermodynamic limit if the thermodynamic potential is a homogeneous 
function of the first order with respect to the extensive variables of state of the system and the entropic 
variable 2 = !/((/— 1) is an extensive variable of state. The equivalence of canonical, microcanonical and 
grand canonical ensembles for the nonrelativistic ideal gas of hadrons was demonstrated. 
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1 Introduction 


The Tsallis power-law distribution describes very well the 
experimental data in both high-energy physics wim and 
astrophysics [1I5] . In particular, the Tsallis-like distribu¬ 
tions are successfully used to describe the transverse mo¬ 
mentum spectra of the particles produced in the proton- 
proton and nucleus-nucleus collisions at the LHC and RHIC 
energies [iEieiniiiniiiiiiiiiis] . However, until now the the¬ 
oretical foundation of the Tsallis statistics mun] is far 
from a definitive solution. The Tsallis generalized statisti¬ 
cal mechanics, by definition, is in the thermal equilibrium. 
Thus, the proof of the thermodynamic self-consistency of 
the Tsallis statistics remains a primordial task, although 
some progress has been made. In this direction, there exist 
some approaches [Biiziiisiiiniiiniiii]. 

The main aim of this paper is to prove analytically 
and numerically the thermodynamic self-consistency of 
the Tsallis statistics in the grand canonical ensemble in 
the approach when the entropic variable 2 : = I/(g — l)is 
considered to be an extensive variable of state and to apply 
this approach of the Tsallis statistics to the hadron ideal 
gas. We investigate the grand canonical ensemble because 
in the high energy collisions the number of particles are 
not exactly conserved and, therefore, the thermodynamic 
quantities and the transverse momentum distributions of 
hadrons in high-energy processes are usually described by 
this statistical ensemble. 


The statistical mechanics defined in a thermal equilib¬ 
rium is thermodynamically self-consistent, i.e., it agrees 
with the requirements of the equilibrium thermodynamics, 
if the thermodynamic potential of the statistical ensemble 
is a homogeneous function of the first order with respect 
to the extensive variables of state of the system in the 
thermodynamic limit The thermodynamic 

self-consistence of the Tsallis statistics in the canonical 
and microcanonical ensembles was proved in pDll2I) . It 
was shown that the homogeneous property of the thermo¬ 
dynamic potential provides the zeroth law of thermody¬ 
namics, the principle of additivity, the Euler theorem and 
the Gibbs-Duhem relation if the entropic variable z is an 
extensive variable of state. Note that the scaling property 
of variable z in the Tsallis statistics and its relation to the 
thermodynamic limit was firstly discussed in |241I25] . 

At present, the concept of the relativistic hadron ideal 
gas is extensively applied in many approaches of high- 
energy physics. It is used to study the hadron production 
mechanism and the QCD phase transition in the heavy-ion 
collisions at high energies However, in the 

present paper we study the thermodynamic properties of 
the hadron gas in the Tsallis statistics in the simple case 
of the hadron gas in the nonrelativistic limit. 

The structure of the paper is as follows. In Section [21 
we formulate the Tsallis statistics in the grand canoni¬ 
cal ensemble. In Sections [3] and 01 we consider the non¬ 
relativistic ideal gas in a finite volume and the themo- 
dynamic limt, respectively. The ensemble equivalence is 
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demonstrated in Section [51 The numerical calculations for 
the hadron ideal gas are given in Section |6l and the main 
conclusions are summarized in the final section. 


2 Tsallis statistics in the grand canonical 
ensemble 

Let us consider the grand canonical ensemble of the dy¬ 
namical systems at the constant temperature T, the vol¬ 
ume T, the chemical potential ^ and the thermodynamic 
coordinate z in a thermal contact with a heat bath. The 
system exchanges with its surroundings the energy and 
matter. The entropy of the Tsallis statistics in the grand 
canonical ensemble is mm 


S = E ^V,ViPiSi, 

(1) 

Si, = kBz{l-pY), 

(2) 


where pi is the probability of the f-th microstate of the 
system, ks is the Boltzmann constant, z = l/{q—l) and 
q GK is a real parameter taking values 0 < q < oo. In 
the limit |z| —>■ oo or q —>• 1, the entropy © recovers the 
Boltzmann-Gibbs entropy, S = —ks J2iPi Note that 
throughout the paper we use the system of natural units 
h = c = ks = i- 

The thermodynamic potential of the grand canonical 
ensemble is related to the fundamental thermodynamic 
potential (energy) by the Legendre transform [50] and can 
be written as 

Q={H)-TS- p{N) = Y, Sv,v,P^n^, (3) 

i 

n, = E,-pN,-Tz(\-pY), (4) 

where Ei and Ni are the energy and number of particles, 
respectively, in the i-th microscopic state of the system. 
The probabilities of microstates in the grand canonical 
ensemble are constrained by an additional function 

‘f = Y ^V’ViPi -1 = 0- (5) 


The unknown probabilities {pi} are obtained from the 
constrained local extrema of the thermodynamic potential 
([3]) by the method of the Lagrange multipliers [3T] 


(p — Q — X(p, 



( 6 ) 

(7) 


where A is an arbitrary real constant. Substituting Eqs. m 
and (I5|) into Eqs. © we obtain 


dHi 

+ Pi -A = 0. 

^Pi 


( 8 ) 


Substituting Eq. (jl]) into ([5]) and using Eq. ®, we ob¬ 
tain the normalized probabilities of the grand canonical 
ensemble as 


Pi = 


1 + 




1 A — Ei -|- pNi 
z + 1 T 

1 A — Ei + pNi 


Vi 


1 + 


z + 1 


T 


(9) 

( 10 ) 


where A = X — T and dEi/dpi = 0. In the limit \z\ —>• 
oo, the probability pi = exp[(yl — Ei + pNi)/T] and A = 
-T\nZ, where Z = J2i exp[-(Ei - pNi)/T]. 

Substituting the probabilities © into Eqs. ©, © and 
using Eq. m we obtain the thermodynamic potential of 
the grand canonical ensemble as 


n = A + T 


1 - 




( 11 ) 


Taking the first derivative of the function f2 with respect 
to temperature T and using Eqs. ®, (nni), and the rela¬ 
tions dEi/dT = 0, dNi/dT = 0 we obtain the entropy 
©-© of the system as 

s =-^ = ( 12 ) 


Applying the Legendre back-transformation to the func¬ 
tion 17 and using Eqs. ©-(HU) and @, we obtain the 
mean energy of the system as 

Bfl 

{H)^n-T— + p{N)^Y^yy^p^^^- (i^) 


Taking the first derivative of the function 17 with respect 
to the chemical potential p and using Eqs. ©, (nni), and 
the relations dEi/dp = 0, dNi/dp = 0 we obtain the 
mean number of particles of the system as 

{N)=-^=Y^VVP^N,- (14) 


The hrst derivative of the function 17 with respect to 
the volume V gives the pressure 


517 


p = 


Ui I, •sr-^ j. 

= -Sv = -Ei''^ 


ViPi 


m 


V dV 


dv 


d5v,Vi 1-1-J 


(15) 


Note that in many applications the last term in Eq. (HU 
is zero. Taking the first derivative of the function 17 with 
respect to the variable z and using Eqs. ©, (nni) and the 
relations dEijdz = 0, dNi/dz = 0 we obtain the conju¬ 
gate force to the variable z as 


A = - 


517 

5z 


TY^yViPi ■ (16) 


% 
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In comparison with the Boltzmann-Gibbs statistics, the 
Tsallis statistics contains a nonvanishing additional con¬ 
jugate force X defined by Eq. (ITCl) pOll^ . 

Equations (HB-dlSl) and the Legendre transform m 
satisfy the differential thermodynamic relation 


df2 = —SdT — pdV — Xdz — {N)dfj. (17) 

and the fundamental equation of thermodynamics mm 

TdS = d{H)+pdV+ Xdz-nd{N). (18) 

The fundamental equation of thermodynamics (1181) pro¬ 
vides the first and second principles of thermodynamics 

5Q = TdS, SQ = d{H) + pdV + Xdz — pd{N), (19) 


where SQ is a heat transfer by the system to the environ¬ 
ment during a quasistatic transition of the system from 
one equilibrium state to a nearby one. Then the heat ca¬ 
pacity can be defined from the general rule, SQ = CdT, 
and the first and second laws of thermodynamics (HU). 
In the grand canonical ensemble the heat capacity at the 
fixed values of (T, V,z,p) can be written as 


Cvzfi 


SQ ds d{H) d(N) 

dT dT dT ^ dT dT^' ^ ^ 


In the grand canonical ensemble the statistical formal¬ 
ism is thermodynamically self-consistent if the thermo¬ 
dynamic potential 12 is a homogeneous function of the 
first degree with respect to the extensive variables of state 
(E,z), i.e., if 


n{T,V,z,p) = Vuj{T,z,,p), (21) 

where Zy = zjV. In this case from Eqs. (fTKI) . (ITfill . (1^ 
and ([3]) we obtain the relation between the pressure and 
the thermodynamic potential and the Euler theorem as 


Substituting Eq. (ESI) into (nni and using the partition 
function of the Maxwell-Boltzmann nonrelativistic ideal 
gas in the Boltzmann-Gibbs statistics 


Z = exp 




27r ) 


we obtain 


Wo ~ jv 
w 


Af=0 


1 + 


1 A^pN' 
JTl T 

3 


z+^N 


"K* =- Ti^} - 


(24) 

, (25) 
(26) 


where oj = gV{mTj2Tr)^^^, g is the spin-isospin degener¬ 
acy factor and m is the particle mass. The upper limit 
Nq in the sum E3 is determined from the conditions 
-z - 3N/2 > 0, 1 -b (A -b pN)/{{z + 1)T) > 0 or/and 
the minimum of the function 


-N 


HN) = —HO) 


1 A + pN' 


-I z+^N 


Z+1 


T 


(27) 


The solution of Eq. EH) is the norm function A as a func¬ 
tion of the variables of state. 

Gombining Eqs. m and (1331) we see that the thermo¬ 
dynamic potential m for the nonrelativistic ideal gas in 
the grand canonical ensemble becomes 


n = A + T 

No - JV 


^ N\' 

N=0 


1 + 


1 yl-b 


z-b 1 T 


z+l+^N 


(28) 


Substituting Eq. (|35)) into the differential equation (ITTl) 
and using Eqs. ED), EH) we have 


uj = —p — Xzy, Ts = e + p + Xzy — pp, (22) 

where uj = fijV, s = S/V, e = {H)/V and p = (N)/V. 
Relations (I3T]) and (1331) can be satisfied in the thermody¬ 
namic limit, i.e. when V oo, \z\ —>• oo and z„ = zjV = 
const. Next, to verify this we shall study a special case of 
this theory, the nonrelativistic ideal gas of hadrons in the 
grand canonical ensemble, in more detail. 


3 Nonrelativistic hadron gas: Grand canonical 
ensemble 

Let us consider the nonrelativistic ideal gas of the Maxwell- 
Boltzmann particles in the framework of the Tsallis statis¬ 
tics in the grand canonical ensemble (T,V, z, p). In this 
paper, we investigate the ideal gas only for the values of 
z < —1. To obtain the exact results, we use the integral 
representation of the Gamma-function |33] 

CXD 

Reix) > 0, Re(y) > 0. (23) 

r{y) J 




w = E 


N=0 


1 + 


1 A + pN 
z-b 1 T 


z+^N 


(29) 


Substituting Eq. (1381) into the differential equation o 
and using Eqs. EH), EHD and EH) we obtain the average 
energy of the system in the grand canonical ensemble as 


No 




N 


A^=0 


(30) 

Making use of Eqs. (1T31) . (1331) . (|33)) and (1351) we can write 
the pressure of the nonrelativistic ideal gas in the grand 
canonical ensemble as 


rp Wo -JV 


V 


N=0 


UJ 

'n\ 


1 A + pN 
z-bl T 


z+l + ^N 


. 

Gomparing Eqs. EQl) and m we obtain the equation of 
state for the nonrelativistic ideal gas in the grand canon¬ 
ical ensemble 

^ 3 V 


0 


(32) 
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Substituting Eq. (EHl) into the differential equation m 
and using Eqs. O we can write the entropy of the system 
in the grand canonical ensemble as 


S = -^{^2-{H)+^^{N)). 


(33) 


As we expected, the sum (12511 in the thermodynamic limit 
dMl) can be approximated by its maximal term, which 
gives the main contribution in (1251) . Thus, from the maxi¬ 
mum of the function (j){N) (l?7l) and Eq. (pK]l we have two 
equations, \n(j){N) = 0 and d\n(j>{N)/dN = 0, which can 
be rewritten as 


This equation coincides with the Legendre transform (|5|). 
Making use of Eqs. (fTHj) and (1^ we can write the conju¬ 
gate force to the variable z in the form 


In ■ 


In^^ 


In^ =0, 


No ~ jV 
LJ' 


N=0 


where 


1 + 


1 A + ^lN 


1 


T 


z+l + ^N 


-In^ 
2 T 


T 


T, 


= 0 


eff 


and 


(34) 


TefJ ^ T 


1 + 


1 Ay + tip 


2 


(40) 

(41) 

(42) 


C = -1 + In 


1 + 


1 A + ^iN' 
z + 1 f 


where w' = Cj/V = g(jnT and A^ = AjV. Here 
(iV) = N corresponds to the maximal term in the sum 
(ED). Erom Eq. (001) we obtain the norm function A as 


- 'ipi-z - 1 - -N) + tp{-z-l). 


(35) 


The mean occupation numbers for the Maxwell-Boltzmann 
nonrelativistic ideal gas in the grand canonical ensemble 
in the framework of the Tsallis statistics can be defined as 


^ n, 

{^ipcr } 


1 + 


1 


pa 


, I 

f'pa • 


n«: 

pa 

1 vi-E 


pa ^pa\z-p 


{Sp p.) 


1 


T 


{iT'pa) — 


we obtain 


A!o ~ Af 
OJ 


{npa) = 


N=0 


1 


1 A + gN - {sp - g) 


z+%N 


z-fl T 

where A is calculated from Eqs. (113 and (123. 


Ay = -fj.p - ZyT 


1-1 + 


2 Zy 


+ // 


T 


Teff = T I 




(43) 

(44) 


where Zq = uj'ejp = {gejp){mT [H]. Substituting 
Eq. (HHl) into Eq. (HT1) and using Eq. obtain the 

equation for the hadron density p as 


(36) 


p = 


— In 
P T 


T, 


eff- 


(45) 


where Sp = p^/2m. Substituting Eq. (123 into (l36ll and 
using the partition function (1241) and the mean occupation 
numbers of the Maxwell-Boltzmann nonrelativistic ideal 
gas in the Boltzmann-Gibbs statistics 


Its solution is the variable p as a function of the variables 
of state {T,Zy,p). Taking Eq. (|25)l in the thermodynamic 
limit (1221) and considering Eq. (HHl) the density of the ther¬ 
modynamic potential can be written as 


(37) 


LU = Ay = —PP — ZyT 


1 - 


1 + ^^ 
2 Zy 


T 


eff 


T 


(46) 


Considering Eqs. (USD, (05]) and (03 we can write the 
energy density in the thermodynamic limit as 


e = -T^ 


d /WN 3 


(47) 


(38) 


Using Eqs. (II3, 03 and (HHl) we obtain the conjugate 
force to the variable z in the thermodynamic limit 


A = - 


duj 

dz. 


= T il- 


T 


eff 


T 


1-ln 


T 


eff 


T 


(48) 


4 Nonrelativistic hadron gas in the 
thermodynamic limit: Grand canonical 
ensemble 

Let us rewrite the thermodynamic quantities of the non¬ 
relativistic ideal gas given above in the thermodynamic 
limit 


V —>■ 00, \z\ —>• 00, Zy 


z 

— = const. 


(39) 


Substituting Eqs. 03, 03 into differential equation 03 
we can write the pressure in the thermodynamic limit as 

duj 2 

p = -u; +Zy-^ = -e = pTeff. (49) 

Making use of Eqs. (1221) . (1221) and (1271) we can write the 
entropy density in the thermodynamic limit as 


uj — £ + pp 

S= - Tf, - = 


1 - 


T 


eff 


T 


(50) 
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Using equation In = 0 we obtain the mean occupa¬ 
tion numbers (l38ll in the thermodynamic limit (1391) as 

(npa) = e . (51) 


5 Ensemble equivalence 

5.1 Equivalence of canonical and grand canonical 
ensembles 


Then, the one-particle distribntion function can be defined 
as 


f{p) 


1 




9 

{2'kYp 



(52) 


where 

J d^Pfip) = 1- (53) 

Let us verify that the nonrelativistic ideal gas in the 
grand canonical ensemble satisfies the principle of addi¬ 
tivity and the zeroth law of thermodynamics in the ther¬ 
modynamic limit. Suppose that the system is divided into 
two subsystems (1 and 2). Then, the extensive variables 
of state of the grand canonical ensemble are additive 


V = V^ + V^, z = z^ + z^ (54) 


and the intensive variables of state of the grand canonical 
ensemble are the same in each subsystem 

T = T^=T\ z,=zl=zl p = p^=p^. (55) 


For the homogeneous functions of the first order, T = 
(T, 17, (iJ), {N),S), and their densities, (p = {Ay,u), e, p, s). 


we have the relation 


F{T,V,z,p) = Vp(T, Zy,p) 

(56) 

and the Euler theorem 


dF dF 

(57) 



Let us demonstrate that the thermodynamic quantities 
of the canonical ensemble in the thermodynamic limit 
can be obtained from the quantities of the grand canon¬ 
ical ensemble by the transformation of the variables of 
state (T,Zy,p,) with the variables of state of the canon¬ 
ical ensemble {T,v,z), where z = z/{N) = Zyjp and 
V = V/(N) = 1/p. Expressing Eqs. (H^ . (H51) in terms 
of the variable z, we obtain the chemical potential of 
the canonical ensemble in the thermodynamic limit (see 
Eq. (38) in pT| l 


9 = 
Teff = T 


5,-1 T'eff 

2+Un — 


T, 


e//’ 


(62) 

(63) 


The energy density (|T7|) of the grand canonical ensem¬ 
ble can be transformed to the energy per particle of the 
canonical ensemble, e = {H)/(N), by using the variable z 
(see Eq. (32) in 


e=- = lTyff. (64) 

p z 

The thermodynamic potential of the canonical ensemble, 
the free energy, can be obtained from the thermodynamic 
potential of the grand canonical ensemble by the Legen¬ 
dre transform, i.e. F = Q + p{N). Then, from Eq. 
we have the thermodynamic potential per particle of the 
canonical ensemble, / = F/{N), which can be written as 
(see Eqs. (30) and (33) in [3J) 


Using Eqs. (l43l) . (l44l) . (1461) . (iTfl) . (l50l) and Eqs. (iMl) . (1^ 
we can verify that the homogeneous functions of the first 
order F are extensive and their densities ip, which are the 
homogeneous functions of the zero order, are intensive 

F = F^+F'^, p = p^=p^. (58) 


/ 


UJ 


-h /i = -zT 

P 




Teff 

T 


(65) 


Rewriting the variable X of the grand canonical en¬ 
semble (I48p in terms of z, we obtain the variable X of the 
canonical ensemble (see Eq. (37) in [21]) 


Eor the homogeneous functions of the zero order, (j) = 
{p,X), we have the relation 


X = T 


1 - 


Teff 

T 


1 - In 



( 66 ) 


c/{T,V,z,p) = c/{T,Zy,p) (59) 


and the Euler theorem 


dV 



= 0 . 


(60) 


Using Eqs. (HSl) . (H51) . and Eqs. (1511) . (1551) we obtain 
that the homogeneous functions of the zero order (/> are 
intensive 


c/ = 4,^ = 4>\ 


(61) 


The pressure of the grand canonical ensemble (091) can be 
rewritten as 


£ = Teff 

3v V 


(67) 


This equation for the pressure of the canonical ensemble 
exactly coincides with Eq. (36) in [3T]. Expressing Eq. (151)1) 
in terms of the variable z we obtain the entropy per parti¬ 
cle of the canonical ensemble, s = S/ (TV), in the form (see 
Eq. (34) in [H]) 


These relations prove the zero law of thermodynamics for 
the Tsallis statistics in the grand canonical ensemble. 


s 


s 

— = Z 


p 



Teff 


( 68 ) 


T 
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Changing the variables of state in Eqs. (ED)-® and 
using Eq. (16211 we obtain the mean occupation numbers 
and the distribution function /(p) for the canonical en¬ 
semble 


(Upfj 



/ mTeff 

V 27r 



(69) 


and 


f{p) = 


1 


2TTmeffT 


3/2 




(70) 


where 



Equations dZOD and (ED exactly coincide with Eqs. (43) 
and (44), respectively, given in [2T] . 

Thus, we have proved that the canonical and grand 
canonical ensembles of the Tsallis statistics are equivalent 
in the thermodynamic limit. 


5.2 Equivalence of microcanonical and grand canonical 
ensembles 


Let us demonstrate that the thermodynamic quantities of 
the microcanonical ensemble in the thermodynamic limit 
can be obtained from the quantities of the grand canonical 
ensemble by the transformation of the variables of state 
(T, Zy, fi) with the variables of state of the microcanonical 
ensemble {s,v,z), where e = {H)/{N) = ejp. Using the 
definition of Zq given in Sec. [Hand Eqs. (H71) . (H^ . we 
can write (see Eq. (48) in [5D]) 



Using again Eqs. (HD and (H^ we obtain 

T = (73) 

Teff = (74) 

Compare Eq. (1751) with Eq. (51) given in Ref. [50]. Sub¬ 
stituting Eqs. dZD, El into Eq. we obtain the en¬ 
tropy per particle for the microcanonical ensemble as (see 
Eq. (49) in [50] ! 


s = 


s 

P 


= z 



(75) 


Considering Eqs. El, El and dH we can write the 
chemical potential for the microcanonical ensemble as (see 
Eq. (54) in Ref. [ID]) 



(76) 


Substituting Eqs. El, dH into Eq. (l4l we obtain (see 
Eq. (55) in Ref. [DD]) 





(77) 


Substituting Eqs. (1751) . (1711) into Eq. (|1D1) we obtain the 
pressure for the microcanonical ensemble as (see Eqs. (53) 
and (51) in Ref. [5D] 1 


T 

p = —w 

V 


- 1 /s 


2 e 

3 V 


(78) 


Thus, we have proved that the microcanonical and 
grand canonical ensembles of the Tsallis statistics are equiv¬ 
alent in the thermodynamic limit. 


6 Analysis and results 

Let us investigate numerically the thermodynamic limit 
for the nonrelativistic ideal gas in the grand canon¬ 
ical ensemble in the framework of the Tsallis statistics. 
We shall increase the volume V at the fixed values of 
the variables of state (T,Zy^fi) and compare the thermo¬ 
dynamic quantities of the ideal gas at constant V with 
their analytical results in the thermodynamic limit. Note 
that here and hereafter po is the normal nuclear density 
and the symbols asterisk and ”G” denote the thermody¬ 
namic quantities in the thermodynamic limit (13911 and in 
the Boltzmann-Gibbs limit, respectively. 

EigurelD shows the dependence of the density of ther¬ 
modynamic potential w, the density of the norm function 
Ay, the entropy density s, the energy density e, the vari¬ 
able X, the reduced particle density p/po and the pressure 
p on the volume V for the nonrelativistic ideal gas of neu¬ 
trons in the grand canonical ensemble in the Tsallis statis¬ 
tics at a given temperature and the chemical potential for 
different values of the entropic variable Zy. For the given 
T, Zy and p, the density of the thermodynamic potential uj 
(I28p and the density of the norm function Ay (HD increase 
with V and attain the constants, ca* = lim cj, ([iS)) 

V^—>-00,2:—)-oo 

and. A* = lim Ay, (155)) . respectively, at very large 

V—)-oo,2—>-oo 

values of V. The density of the thermodynamic potential 
oj is much lower than the density of the norm function 
Ay at finite values of V and it tends to Ay with V. In 
the thermodynamic limit as U —oo and at Zy constant 
the function Ay achieves the finite value w*; however, the 
density of the thermodynamic potential cu* differs essen¬ 
tially from the density of the thermodynamic potential of 
the Boltzmann-Gibbs statistics, = lim uj*. At finite 

2.L,—>-00 

values of and large values of the volume V the den¬ 
sity of the norm function A* coincides with the density 
of the thermodynamic potential uj* and differs from its 
Boltzmann-Gibbs limit uj^. The density of the thermody¬ 
namic potential oj* increases with \zy\ at fixed variables 
of state (T, p) and attains its Boltzmann-Gibbs value uj^ 
at very large values of \zy \ —>■ oo. The functions oj and Ay 
tend to minus infinity with decreasing the volume V. 

The shape of the entropy density in the thermody¬ 
namic limit is shown in the middle left panel of Fig. [D The 
entropy density s (1551) decreases with V at fixed values of 
(T, Zy, p) and attains the constant, s* = lim s, (I5D1) 

V —>-00,2—foo 

as U —>■ 00 . At finite values of Zy, the entropy density 
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Fig. 1. (Color online) The density of the thermodynamic potential u), the density of the norm function yl„, the entropy density 
s, the energy density e, the variable X, the density p and the pressure p as functions of the volume V for the nonrelativistic 
ideal gas of neutrons in the grand canonical ensemble in the Tsallis statistics at temperature T = 100 MeV and the chemical 
potential p = —200 MeV for different values of The curves 1, 2, 3 and 4 correspond to the values of Zv = —5, —10, —20 and 
— 100 fm“^, respectively. The dashed lines represent the Tsallis statistics in the thermodynamic limit. The dotted lines represent 
the Boltzmann-Gibbs statistics and the dash-dotted lines correspond to the quantity /l„. 


in the thermodynamic limit s* differs essentially from the 
entropy density of the Boltzmann-Gibbs statistics = 
lim s*. The entropy density s* decreases with \zy\ and 

Zv—^OO 

achieves the Boltzmann-Gibbs limit only as \zy \ —1 oo. 

The behaviour of the energy density in the thermody¬ 
namic limit is shown in the bottom left panel of Fig. [TJ 
The energy density s dSO]) decreases with V at fixed val¬ 
ues of the variables of state (T, Zy, p) and attains the con¬ 
stant value, £* = lim e, (SZl) at y 00 . At finite 
V —¥CCi^Z—¥00 

values of the energy density in the thermodynamic 
limit £* differs from the energy density of the Boltzmann- 
Gibbs statistics = lim e*. The energy density e* de- 

Z.u—¥00 

creases with \zy\ and achieves the Boltzmann-Gibbs limit 
as \zy \ —>■ 00 . 


The variable X ((Mil as a function of the volume V at 
fixed values of the variables of state (T, Zy, p) is shown in 
the upper right panel of Fig. [T] It also decreases with V 
and attains the constant value, X* = lim X, (H51) 

V —¥oo,z—>co 

at y —>■ 00 . For the finite values of Zy, the function X* is 
nonvanishing and differs from its Boltzmann-Gibbs zero 
value, X^ = Q. The conjugate force X* decreases with 
\zy \ and vanishes as |z„| —oo. 

The density of neutrons (ESI) is shown in the middle 
right panel of Fig. [T] It decreases with V and attains the 
constant value p* (1151) as y ^ oo. For finite values of Zy, 
the neutron density p* differs from its Boltzmann-Gibbs 
limit p^. The density p* decreases with \zy\ and achieves 
the Boltzmann-Gibbs limit p'^ as \zy \ —> oo. 

The shape of the pressure (15T1) is shown in the bottom 
right panel of Fig. [T] It decreases with V and attains the 
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Fig. 2. (Color online) The shift = uj* — to and the shift 
— oj as functions of the volume V for the non- 
relativistic ideal gas of neutrons in the grand canonical en¬ 
semble of the Tsallis statistics at the temperature T = 100 
MeV and chemical potential /r = —200 MeV for different val¬ 
ues of Zv The curves 1, 2, 3 and 4 correspond to the values of 
Zv = —5, —10, —20 and —100 fm“®, respectively. 

constant value p* (PI) as y —>■ 00 . For finite values of 
Zy, the pressure p* differs from its Boltzmann-Gibbs limit 
p^. The pressure p* decreases with \zy\ and achieves the 
Boltzmann-Gibbs limit p^ as \zy \ —>■ oo. 

Figure [5] illustrates the dependence of the shift = 
oj* — oj and Ai^a = Ay — cu on the volume V for the non- 
relativistic ideal gas of neutrons in the grand canonical 
ensemble in the Tsallis statistics at the given temperature 
and chemical potential for different values of the entropic 
variable Zy. At the fixed values of the variables of state 
{T,Zy,pL), the deviation of the density of the thermody¬ 
namic potential u from its thermodynamic limit value w* 
decreases with V and it vanishes as —>■ oo. The be¬ 
haviour of the shift A^j with V confirms numerically the 
thermodynamic limit (j39ll for the Tsallis statistics in the 
grand canonical ensemble. The shift A^^a also decreases 
with V and vanishes as V oo. Such behavior proves 
numerically the equality of the density of the thermody¬ 
namic potential w* with the density of the norm function 
Ay in the thermodynamic limit. 

Figure [3] shows the dependence of the shift Ap = p — 
p* and the shift A^ = e — e* on the volume V for the 
nonrelativistic ideal gas of neutrons in the grand canonical 
ensemble of the Tsallis statistics at the given temperature 
T and chemical potential p for different values of Zy. The 
functions Ap and Ag decrease with V and tend to zero as 
V ^ oo. Such behaviour of the shifts Ap and A^ proves 



10 ° 10 ' 10 ° 10 ° 10 '' 10 ° 
V(fTn°) 


Fig. 3. (Color online) The shift Ap = p — p* and the shift 
Ae = e — e* as functions of the volume V for the nonrela¬ 
tivistic ideal gas of neutrons in the grand canonical ensem¬ 
ble of the Tsallis statistics at the temperature T = 100 MeV 
and chemical potential p = —200 MeV for different values 
of Zy. The curves 1,2,3 and 4 correspond to the values of 
Zy = —5, —10, —20 and —100 fm~®, respectively. 

numerically the thermodynamic limit (1391) for the density 
p and the energy density e of the Tsallis statistics. 

Let us verify numerically the homogeneity properties 
of the thermodynamic potential of the nonrelativistic ideal 
gas of neutrons in the grand canonical ensemble for the 
Tsallis statistics. We shall find the range of the variables 
of state V and z at the fixed values of (T, Zy , p) where 
the thermodynamic quantities of the Tsallis statistics sat¬ 
isfy the main requests of the equilibrium thermodynamics. 
The formalism of the statistical mechanics in the grand 
canonical ensemble agrees with the requirements of the 
equilibrium thermodynamics, i.e. it is thermodynamically 
self-consistent, if the thermodynamic potential of the grand 
canonical ensemble is a homogeneous function of the first 
order with respect to the extensive variables of state. Thus, 
the Tsallis statistics in the grand canonical ensemble will 
be thermodynamically self-consistent if its thermodynamic 
potential is a homogeneous function of the first order with 
respect to the extensive variables of state V and z. This 
means that the density of the thermodynamic potential w 
for the grand canonical ensemble shall satisfy the following 
equation; 

A = uj p Zy = 0 

for any values of the variables of state (T, Zy, p). This equa¬ 
tion is a result of Eqs. m and (1221) . 

Let us verify numerically Eq. (1731) for the ideal gas of 
the Tsallis statistics in a finite volume. Figure 0] displays 
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Fig. 4. (Color online) The shift A = ui -\- p -\- Xzv as a func¬ 
tion of the volume V for the nonrelativistic ideal gas of neu¬ 
trons in the grand canonical ensemble of the Tsallis statis¬ 
tics at the temperature T = 100 MeV and chemical potential 
p = —200 MeV for different values of z„. The curves 1,2,3 
and 4 correspond to the values of z-u = —5,—10,—20 and 
— 100 fm“®, respectively. The dashed line represents the shift 
A* = uj* + p* + X*Zv in the thermodynamic limit. 


the dependence of the shift A on the volume V for the 
nonrelativistic ideal gas of neutrons in the grand canonical 
ensemble of the Tsallis statistics at the given temperature 
T and chemical potential /i for different values of It 
is clearly seen that at small values of the volume V the 
shift A is not equal to zero and the density of the ther¬ 
modynamic potential uj is an inhomogeneous function. At 
large values of volume V in the thermodynamic limit the 
shift A vanishes and the density of the thermodynamic 
potential lo becomes a homogeneous function of the first 
order with respect to the extensive variables of state V 
and 2 . Thus, we have demonstrated numerically that the 
Tsallis statistics in the grand canonical ensemble becomes 
thermodynamically self-consistent in the thermodynamic 
limit (I55|) . 

Figure [5] represents the dependence of the one-particle 
distribution function (13811 on the momentum p of a particle 
for the nonrelativistic ideal gas of neutrons in the grand 
canonical ensemble of the Tsallis statistics at the given 
temperature T, chemical potential p, and variable for 
different values of V and z. At small values of the vol¬ 
ume V and the variable |z„|, where the Tsallis statistics is 
thermodynamically inconsistent, the one-particle distribu¬ 
tion function of the Tsallis statistics has a power-law form 
and its shape deviates essentially from the usual Maxwell- 



p (GeV) 

Fig. 5. (Color online) The dependence of the one-particle dis¬ 
tribution function on the momentum p for the nonrelativistic 
ideal gas of neutrons in the grand canonical ensemble of the 
Tsallis statistics at the temperature T = 100 MeV and chemi¬ 
cal potential p — —200 MeV for two values of Zv and different 
values of V and 2 . The solid curves 1, 2 and 3 correspond to 
the values of V = 1.728,3.375 fm® and V = oo, respectively. 
The dashed lines represent the Maxwell-Boltzmann distribu¬ 
tion function of the Boltzmann-Gibbs statistics. 


Boltzmann distribution function. The Tsallis one-particle 
distribution function /(p) tends to the distribution func¬ 
tion /*(p) in the thermodynamic limit as V —1 oo. The 
Tsallis distribution function in the thermodynamic limit 
/*(p) has an exponential form and it is close to the Maxwell- 
Boltzmann distribution function of the Boltzmann-Gibbs 
statistics, but they do not coincide. The difference be¬ 
tween the distribution functions /(p) and /*(p) dimin¬ 
ishes with \zv\. For small values of the variable \zv\, the 
Tsallis one-particle distribution function in the thermody¬ 
namic limit f*{p) can differ essentially from the Maxwell- 
Boltzmann distribution function f'^{p) of the Boltzmann- 
Gibbs statistics. 

Figure [5] shows the dependence of the density of the 
thermodynamic potential w*, the entropy density s*, the 
energy density e* , the variable X* , the density p* and the 
pressure p* for the Tsallis statistics in the thermodynamic 
limit on the variable \zy\ for the nonrelativistic ideal gas 
of neutrons in the grand canonical ensemble at the given 
chemical potential p for different values of temperature T. 
At small values of |z«|, the thermodynamic quantities w*, 
s*, £*, X*, p* and p* of the Tsallis statistics in the ther¬ 
modynamic limit differ essentially from their values in the 
Boltzmann-Gibbs limit. Otherwise, at large values of \zy\, 
as \zv\ oo, these quantities resemble the corresponding 
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Fig. 6. (Color online) The Boltzmann-Gibbs limit for the Tsallis statistics in the thermodynamic limit. The density of the 
thermodynamic potential uj*, the entropy density s*, the energy density e*, the variable X*, the density p* and the pressure p* 
for the Tsallis statistics in the thermodynamic limit as functions of the variable Zy for the nonrelativistic ideal gas of neutrons in 
the grand canonical ensemble at the chemical potential p = —200 MeV for different values of temperature T. The curves 1 and 
2 correspond to the values of T = 100 and 110 MeV, respectively. The dashed lines represent the Boltzmann-Gibbs statistics. 


thermodynamic quantities of the Boltzmann-Gibbs statis¬ 
tics. The entropy density, the energy density, the neutron 
density and the pressure for both the Tsallis statistics in 
the thermodynamic limit and the Boltzmann-Gibbs statis¬ 
tics increase with temperature T. However, the density of 
the thermodynamic potential for both the Tsallis statis¬ 
tics in the thermodynamic limit and the Boltzmann-Gibbs 
statistics decreases with T. 

Figure [7] represents the dependence of the density of 
the thermodynamic potential lo*, the entropy density s*, 
the energy density £*, the variable X*, the density p* and 
the pressure p* for the Tsallis statistics in the thermo¬ 
dynamic limit on the temperature T for the nonrelativis¬ 
tic ideal gas of neutrons in the grand canonical ensem¬ 
ble at the given chemical potential p for different val¬ 
ues of the variable At small values of the tempera¬ 
ture T and the fixed values of |z„|, the Tsallis statistics 


in the thermodynamic limit approaches the Boltzmann- 
Gibbs statistics. However, at large values of the tempera¬ 
ture T, the thermodynamic quantities of the Tsallis statis¬ 
tics in the thermodynamic limit differ from their values 
in the Boltzmann-Gibbs limit and this deviation increases 
with T. For the negative values of the variable Zy, the ther¬ 
modynamic quantities s*, e*, X*, p* and p* of the Tsal¬ 
lis statistics in the thermodynamic limit overestimate the 
corresponding quantities of the Boltzmann-Gibbs statis¬ 
tics; however, the density of the thermodynamic potential 
underestimates the density of the thermodynamic poten¬ 
tial of the Boltzmann-Gibbs statistics. 

Figure [5] shows the dependence of the density of the 
thermodynamic potential w*, the entropy density s*, the 
energy density £*, the variable X*, the density p* and 
the pressure p* for the Tsallis statistics in the thermo¬ 
dynamic limit on the chemical potential p for the non- 
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Fig. 7. (Color online) Comparison of the Boltzmann-Gibbs statistics with the Tsallis statistics in the thermodynamic limit. 
The density of the thermodynamic potential uj*, the entropy density s*, the energy density e*, the variable X*, the density p* 
and the pressure p* for the Tsallis statistics in the thermodynamic limit as fnnctions of the temperature T for the nonrelativistic 
ideal gas of neutrons in the grand canonical ensemble at the chemical potential p = —200 MeV for different values of the variable 
Zv The solid curves 1 and 2 correspond to the values of Zv = —5 and —10 fm“®, respectively. The dashed curves represent the 
Boltzmann-Gibbs statistics. 


relativistic ideal gas of neutrons in the grand canonical 
ensemble at the given temperature T for different values 
of the variable Zy. For the fixed values of \zy\, the Tsal¬ 
lis statistics in the thermodynamic limit approaches the 
Boltzmann-Gibbs statistics at large absolute values of the 
chemical potential p. However, at small absolute values of 
the chemical potential p, the thermodynamic quantities 
of the Tsallis statistics in the thermodynamic limit differ 
from their values in the Boltzmann-Gibbs limit and this 
deviation increases with |/i| 0. For the negative values 

of the variable the thermodynamic quantities s*, e*, 
X*, p* and p* of the Tsallis statistics in the thermody¬ 
namic limit overestimate the corresponding quantities of 
the Boltzmann-Gibbs statistics; however, the density of 
the thermodynamic potential underestimates the density 
of the thermodynamic potential of the Boltzmann-Gibbs 
statistics. 


7 Conclusions 

In the present paper, the Tsallis statistics in the grand 
canonical ensemble was formulated in the general form. 
The thermodynamic potential of the grand canonical en¬ 
semble was found from the fundamental thermodynamic 
potential by the Legendre transform. The probabilities 
of microstates of the system were derived from the con¬ 
strained local extrema of the statistical thermodynamic 
potential of the grand canonical ensemble by the method 
of the Lagrange multipliers. It was proved that both inten¬ 
sive and extensive thermodynamic quantities of the Tsallis 
statistics satisfy the same differential thermodynamic re¬ 
lations and statistical definitions in terms of probabilities 
Pi as those for the Boltzmann-Gibbs statistics. This is a 
consequence of the thermal equilibrium and the Legendre 
transform. It was shown that in the Tsallis statistics the 
extensive variable of state 2 generates the nonvanishing 
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Fig. 8. (Color online) Comparison of the Boltzmann-Gibbs statistics with the Tsallis statistics in the thermodynamic limit. 
The density of the thermodynamic potential uj*, the entropy density s*, the energy density e*, the variable X*, the density 
p* and the pressure p* for the Tsallis statistics in the thermodynamic limit as fnnctions of the chemical potential p for the 
nonrelativistic ideal gas of neutrons in the grand canonical ensemble at the temperature T = 100 MeV for different values of 
the variable Zv The solid curves 1 and 2 correspond to the values of Zv = —5 and —10 fm“®, respectively. The dashed curves 
represent the Boltzmann-Gibbs statistics. 


conjugate force X. Moreover, the fundamental equation of 
thermodynamics was consistently derived from the ther¬ 
modynamic potential of the grand canonical ensemble. 

The exact analytical relations for the thermodynamic 
quantities (thermodynamic potential and its first deriva¬ 
tives) of the nonrelativistic ideal gas in the framework of 
the Tsallis statistics in the grand canonical ensemble were 
obtained. They were analytically rewritten in the thermo¬ 
dynamic limit considering the variable z as an extensive 
variable of state. Based on the example of the nonrela¬ 
tivistic ideal gas in the grand canonical ensemble it was 
proved that the Tsallis statistics satisfies the requirements 
of the equilibrium thermodynamics in the thermodynamic 
limit when the entropic variable z is an extensive variable 
of state of the system, i.e., the thermodynamic potential is 
a homogeneous function of the first order with respect to 
the extensive variables of state of the system. The equiva¬ 


lence of the canonical, grand canonical and microcanoni- 
cal ensembles in the thermodynamic limit was analytically 
proved. It was shown that in the thermodynamic limit the 
relations of the ideal gas of the Tsallis statistics in terms 
of the function can be written in a form similar to the 
Boltzmann-Gibbs statistics. The behaviour of the thermo¬ 
dynamic quantities of the ideal gas in both the thermo¬ 
dynamic limit and the Boltzmann-Gibbs limit were also 
studied. It was numerically proved that the thermody¬ 
namic potential of the grand canonical ensemble becomes 
a homogeneous function of the first order with respect to 
the extensive variables of state only in the thermodynamic 
limit. It was shown that the one-particle distribution func¬ 
tion has a power-law form when the Tsallis statistics is 
thermodynamically inconsistent and takes an exponential 
form in the thermodynamic limit when the Tsallis statis¬ 
tics is thermodynamically consistent. 
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